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ABSTRACT 

We present mid-infrared (5.2-15.2 fim) spectra of the Type la supernovae (SNe la) 2003hv and 
2005df observed with the Spitzer Space Telescope. These are the first observed mid-infrared spectra 
of thermonuclear supernovae, and show strong emission from fine-structure lines of Ni, Co, S, and Ar. 
The detection of Ni emission in SN 2005df 135 days after the explosion provides direct observational 
evidence of high-density nuclear burning forming a significant amount of stable Ni in a Type la 
supernova. The observed emission line profiles in the SN 2005df spectrum indicate a chemically 
stratified ejecta structure. The Ar and Ni lines in SN 2005df are both shifted to the red relative 
to the Co emission, which appears to be nearly centered with respect to the rest frame of the host 
galaxy. The SN 2005df Ar lines also exhibit a two-pronged emission profile implying that the Ar 
emission deviates significantly from spherical symmetry. The observed Ar lines can be reproduced 
with either an almost circular ring-shaped geometry observed nearly edge-on, or a prolate ellipsoidal 
emission geometry with a central hole observed nearly pole-on. The spectrum of SN 2003hv also 
shows signs of asymmetry, exhibiting blueshifted [Co III] which matches the blueshift of [Fe II] lines 
in nearly coeval NIR spectra. Finally, local thermodynamic equilibrium abundance estimates for the 
yield of radioactive ^^Ni give Mse^i ~ 0.5Mq, for SN 2003hv, but only Msejv, « O.13-O.22M0 for the 
apparently subluminous SN 2005df, supporting the notion that the luminosity of SNe la is primarily 
a function of the radioactive ^^Ni yield. 

The chemically stratified ejecta structure observed in SN 2005df matches the predictions of delayed- 
detonation (DD) models, but is entirely incompatible with current three-dimensional deflagration 
models. Furthermore the degree that this layering persists to the innermost regions of the supernova 
is difficult to explain even in a DD scenario, where the innermost ejecta are still the product of 
deflagration burning. Thus, while these results are roughly consistent with a delayed detonation, it is 
clear that a key piece of physics is still missing from our understanding of the earliest phases of SN la 
explosions. 

Subject headings: supernovae: general — supernovae: individual (SN 2003hv; SN 2005df) 



1. INTRODUCTION 

There is general agreement about the basic nature of 
Type la supernovae (SNe la): they are the thermonu- 
clear explosion of degenerate white dwarfs (WDs), prob- 
ably in some sort of close binary system undergoing 
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some form of mass transfer, and exploding once a criti- 
cal mass is reached, likely near the Chandrasekhar limit, 
at least for the majority of SN la events. (For reviews , 
see Branch 1998; Lcibundgut 2001; Nom oto et 311120031 : 
IHoflich et al.1[2003l ) StiU, many of the details of this pic- 
ture remain either unknown or are highly controversial. 
In particular, the initial conditions of the WD progeni- 
tor prior to thermonuclear runaway, and the dynamical 
propagation of the nuclear burning are both hotly de- 
bated topics of discussion. 

The recent growth in the availability of infrared (IR) 
observations shows much promise in addressing these 
issues. Infrared spectroscopy has several fundamental 
advantages over UV/optical spectroscopy for work on 
SNe la. First, the strong spectral features in the IR are 
spread farther apart in wavelength, avoiding most of the 
strong line-blending issues that complicate the analysis of 
optical/UV observations. Consequently, kinematic infor- 
mation can be extracted directly from the observed line 
profiles without relying on a detailed spectral synthesis 
calculation. Second, IR lines tend to have a lower op- 
tical depth than UV/optical lines which, combined with 
the reduced line blending, greatly aids in the accurate 
measurement of abundances of iron-peak elements. 

This is particularly true of ground-state transitions in 
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the mid-infrared (MIR). For these low-excitation lines 
the effects of stimulated emission are very strong, dra- 
matically reducing the optical depth. Thus, these lines 
can become optically thin quite rapidly, often while the 
electron densities are still well above the critical density 
for coUisional de-excitation. As a result, relatively accu- 
rate (and temperature insensitive) abundance estimates 
can be made from MIR lines using just a local thermody- 
namic equilibrium (LTE) approximation. Finally, near- 
IR and mid-IR observations can be used to put strong 
constraints on the formation of molecules and dust in the 
cooling ejecta. 

However, IR observations have one disadvantage. 
While the IR flux of SNe la is intrinsically fainter than 
in the optical, the foreground and background emission 
is usually much larger. In particular, foreground thermal 
emission precludes MIR spectroscopy of all but the very 
nearest extragalactic SNe la from the ground, and to 
date, IR spectroscopic observations of SNe la have been 
limited to the near-infrared ( NIR). With the lau nch of 
the Spitzer Space Telescope (| Werner et al.ll2004) . MIR 
spectroscopy is now possible, at least for nearby {D < 
25 Mpc) SNe la. In this paper we present the first 
MIR spectra of SNe la, obtained w ith the Infrared Spec- 
trograph fIRS: [Houck et al]|2004D on the Spitzer Space 
Telescope: A low signal-to- noise ratio (S/N) spectrum of 
SN 2003hv and the first of two planned high S/N obser- 
vations of SN 2005df. Both were obtained as part of the 
GO program of the Mid-Infrared Supernova Gonsortium 
(MISC). 

In § 2 we describe the observations and data reduction, 
and the results are shown in § 3, including a detailed dis- 
cussion of line identifications. In § 4 we analyze the ob- 
served emission-line profiles using synthetic line profiles 
from simple models of the emission distribution. In § 5, 
we use the observed nickel and cobalt line fluxes to infer 
abundances assuming LTE level populations, and make 
estimates for the optical depth and electron density. Fi- 
nally, we discuss and summarize our results in § 6 and 
§7. 

2. OBSERVATIONS AND DATA REDUCTION 

SN 2003hv was discovered at about 12. 5 mag on 9.5 
Sep 2 003 (all dates UT) in NGG 1201 by iBeutler fc Lil 
(120031). and cl assified as a SN la on 10.3 Sep. 2003 by 
iDressler et all ^2003), closely resembhng SN 1994D two 
days after maximum light. The Lyon/Meudon Extra- 
galactic Database (LEDA)^'^ gives a redshift corrected 
for Virgo infall of 1494 km s~^. Assuming a Hubble con- 
stant iJo = 70 km s~^ Mpc~^, this implies a distance of 
21.3 Mpc. No light-curve data have yet been published 
for SN 2003hv, but the reported discovery magnitude , 
Galactic extinction of Ay = 0.05 (|Schlegel et al.lll998f ). 
and a 21.3 Mpc distance implies a peak magnitude of 
My ~ —19.2. Together with the reported initial spec- 
trum, this suggests that SN 2003hv was a relatively nor- 
mal Type la event. Assuming that Bmax is indeed two 
days earlier than the initial spectrum (8.5 Sep. 2003), 
and SN 2003hv is a no rmal SN la, then a 19.5 day rise 
time (Riess "et~aIIIT999l) gives an estimated explosion date 
of 22 Aug. 2003. (Note that throughout this paper we 
will date the SN observation epochs with < = at the 



explosion, rather than the more common date of i?max-) 
SN 2005df was discove r ed vis ually in NGG 1559 on 
4 Aug. 2005 by lEvand ()2005D . and co nfirmed on 5 
Aug. 2005 with a CCD image bv iGilmorc ( 20051) It 
was c lassified as a peculiar SN la bv iSalvofc Schmidtl 
()2005[ ). who reported that its optical spectrum on 5 
Aug. 2005 resembled that of the moderately sublumi- 
nous Type la SN 1986G a few days before maximum 
light. As with SN 2003hv, no light-curve data have yet 
been published, but observations reported on SNWeb^** 
suggest that SN 2005df peaked near 12.3 mag around 18 
Aug. 2005. 

Several different distance estimates to NGC 1559 have 
been given. The LEDA redshift of NGG 1559 corrected 
for Virgo infall is 1005 km s~^, which yields a dis- 
tance of 14.4 Mpc for Hq = 70 km s"^ Mpc"^ An 
expanding photosph e re method analysis of SN 1986L 
(jSchmidt et all 11994 lEastman et"aI1 Il996l) gives a dis- 
tance of 16 ±2 Mpc, and lHamuvl (l2002l) gives an e stimate 
of 18.7 ±2.2 Mpc based on the lTonrv et all ()2000f ) model 
of the peculiar velocity flow. Adopting a p eak magnitude 
of 12 .3, Galactic extinction Ay = 0.1 mag (jSchlegel et al.l 
Il998| ). and no host-galaxy extinction yields peak lumi- 
nosities of My « -18.5 and -19.0 mag for D = 14.4 Mpc 
and 18.7 Mpc, respectively. This suggests that SN 2005df 
was subluminous by ~ 0.5-1 mag, relative to nor mal 
SNe la with My « -19.5 mag (|Gibson et al.ll2000l ). It 
is perhaps possible that some of this is due to extinction 
in the host galaxy, but F-band and i?-band observations 
reported on SNWeb give {V — R) w mag, suggesting 
that SN 2005df is subject to relatively little reddening. 
Assuming Vmax occurred on 18 Aug., and adopting a 16 
day rise time from delayed-detonation mo dels appropri- 
ate t o a moderately subluminous SN la (jHoflich et al.l 
[200l) . gives an estimated explosion date of 2 Aug. 2005. 

Spitzer observations of SN 2003hv consisted of a sin- 
gle set of short-low (5.2-15.2 /xm) IRS spectra, obtained 
as part of the GOl program of the MISC. Observations 
of SN 2005df were obtained via a low-impact target-of- 
opportunity (ToO) trigger as part of the MISG G02 
program. The full set of ToO observations includes 
two epo chs of imag i ng wi th the Infrared Array Camera 
(IRAC; iFazio et al.ll2004l ). two epochs of imaging with 
the blue peak- up array of the IRS (PUI), and 1 epoch 
of short-low spectroscopy with IRS. In addition to these 
ToO observations, a second epoch of both short-low and 
long-low (14-38 /im) IRS spectroscopy and a third set of 
IRAC/PUI imaging is scheduled for observation in G03. 
The full set of Spitzer data on SN 2005df will be analyzed 
in a later paper. Here we will focus primarily on the first 
epoch of IRS spectroscopy. 

IRAC and PUI images of SN 2005df were obtained at 
three epochs, in Nov. 2005, Feb. 2006, and Aug. 2006, 
with the first two bracketing the IRS spectral observa- 
tions. The imaging observations usually consisted of 30 s 
exposures at 20 dither positions for each band, except for 
the Aug. 2006 PUI images, for which four 30 s exposures 
were obtained for each of the 20 dither positions. 

Aperture photometry was performed on the Post-Basic 
Calibrated Data (PBGD) mosaic images using standard 
IRAF^^ routines and a circular aperture with a three- 
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pixel radius centered on the measured centroid of the 
detected point source image. The background was es- 
timated using a circular annulus with inner and outer 
radii of three and seven pixels, respectively. The result- 
ing fluxes were multiplied by the aperture correction fac- 
tors listed in version 3.0 of the IRAC Data Handbook 
and version 2.0 of the IRS Data Handbook. Photometric 
errors due to both photon statistics and background un- 
certainties (estimated using the standard deviation in the 
background annulus) were calculated, with the latter by 
far the dominant source of uncertainty. For images with- 
out a clear detection of the SN, an upper limit to the flux 
was estimated using the uncertainty in the background 
emission at the location of the SN. 

IRS spectroscopy took place on 1 Sep. 2004 and 14 
Dec. 2005 for SN 2003hv and SN 2005df, respectively 
Light curves for these events have not yet been published, 
so we will adopt our above estimates for the explosion 
dates placing the epochs of IRS observation at 375 d for 
SN 2003hv and 135 d for SN 2005df. These epochs are 
probably accurate to a few days, and this degree of timing 
uncertainty does not significantly affect our analysis. 

IRS observations consisted of four sets of exposures: 
two nod positions spaced 19" apart in each of the two 
short-low slits (nominally 1st and 2nd orders). Fifteen 
60 s exposures were obtained at each nod position in each 
order, resulting in a total of 30 min integration time for 
each spectral order. 

Spectral reductions began with the Basic Calibrated 
Data (BCD) products from the SSC pipeline (version 
S13.0.1). The fifteen exposures for each nod position 
were combined using standard IRAF routines. Back- 
ground sky subtraction was then achieved by subtract- 
ing the first nod position in each slit from the second. 
These subtracted frames were then cleaned using the con- 
tributed IRSCLEAN_MASK software package^'^, to re- 
move hot pixels near the SN spectrum. One-dimensional 
(1-D) flux-calibrated spectra were then extracted from 
the cleaned 2-D frames with the S13 version of SPICE^^ 
using the mask output from IRSCLEAN_MASK and the 
uncertainty frames from the PBCD data products. 

The output of this process was four sets of flux- 
calibrated 1-D spectra, one from each nod position in 
each slit. In addition, each of the extracted second- 
order spectra included data from the first-order "bonus 
segment" at the end of the detector. So in total, the 
extracted short-low data consisted of 6 sets of partially 
overlapping 1-D spectra, two each covering 5.2-7.6 /im 
(second order), 7.3-8.7 ^m ("bonus segment") and 7.5- 
15.2 /j,m (first order). 

These six 1-D spectra were then combined using cus- 
tom software to make a variance- weighted average of the 
data from each overlapping data set on a bin-by-bin ba- 
sis, using the error spectrum output by SPICE, and elim- 
inating any data points with abnormal flags. In the over- 
lap region, data from the different spectral orders were 
combined by including flux from partially overlapping 
bins in proportion to the size of the overlap compared 
with the size of the bin (essentially assuming that the 

servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^^^http: //ssc . spitzer.caltech.edu/archanaly/contrlbuted/lrsc 
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flux is spread evenly across a given spectral bin). This 
procedure was found to significantly increase the S /N of 
the resulting combined spectrum compared with just a 
straight average of the overlapping orders. This was par- 
ticularly true in the region near 7 /im, where it is nec- 
essary to combine data near the ends of the individual 
spectral orders because the S/N declines rapidly. 

3. RESULTS 

3.1. SN 2005df Photometry 

The broad-band photometric fluxes from the MIR im- 
ages are presented in Tableland are plotted in Figure [TJ 
The fluxes have relatively large uncertainties, particu- 
larly at the longer wavelengths, and these are, by far, 
dominated by the fluctuations in the background galaxy 
emission near the SN. In a subsequent paper we will try 
to address this with a more careful extraction of the SN 
fluxes, but for now we will restrict ourselves to a fairly 
qualitative discussion of the photometry. 

SN 2005df fades monotonically in all the observed 
bands, but for the most part there is little evidence of 
significant evolution in the observed spectral energy dis- 
tribution (SED). The one clear exception to this is the 
rapid decline in the 3.6 fim (IRAC Channel 1) flux, which 
is high in the first epoch, but has fades by nearly a factor 
of 4 in the 100 days between the first two epochs. 

The large 3.6 /im/4.5 /iin ratio in the first epoch sug- 
gests that strong line emission rather than continuum 
dominates the 3.6 /im flux. Potential candidate emis- 
sion features in this band include [Fe III] 3.229 /im, and 
[Co III] 3.492 /im, though L-band spectra of a SN la 
around 100 days are needed to make a secure identifica- 
tion. Both of these lines have somewhat higher excitation 
temperatures than the spectral features likely dominat- 
ing the 8.0 /im and 16 /im bands, and thus the more rapid 
fading of the 3.6 fim band could be due to cooling in the 
ejecta. A large number of [Fe II] and [Co II] lines also 
lie in the 3.6 /im band but, as we will discuss in S 5.2, 
it seems likely that doubly ionized iron-peak species are 
dominant in the ejecta at these epochs. 

The flux is quite low at 4.5 /im and 5.8 /im at all 
epochs, the latter consistent with the lack of emission 
seen at the blue end of the IRS spectrum (Fig.[2l § 3.2). 
This is in contrast to SNe II, which often show a strong 
flux excess at these wavelengths due to emission from 
the f undamental ro- vibratio nal band of carbon monoxide 
(CO; iCatchpole etlll 119881: Woodc n|[l98l iKotak et all 
|200l 120061). The M-band ( Catchpole et al.lll988f ) fliK 
from SN 1987A at 200 d (roughly the peak of CO emis- 
sion) would correspond to about 1.5-2 mJy at the dis- 
tance (14.4-18.7 Mpc) of SN 2005df. All things being 
equal, the 4.5 /im fluxes of SN 2005df would thus im- 
ply an upper limit on the CO formation of < 5% of that 
formed in SN 1987A. In practice, a more detailed analysis 
is needed to account for radiation transport effects, and 
for the very different physical conditions in the ejecta of 
SNe la and SNe II. 

SN 2005df also shows significant emission at 8.0 fim 
and 16 /im. Our IRS spectrum suggests that the 8.0 /im 
emission is dominated by strong [Ar II] and [Ar III] emis- 
sion, with some contribution from nickel lines (see § 3.2). 
Possible lines contributing to the 16 /im emission include 
le^p II] 14.740, 15.647 and 16.30 /im [Co III] 16.391 /im, 
[Co IV] 15.647 fim (the ground state fine-structure line 
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of Co IV) and [Fe II] 17.936 fiui. Subsequent IRS "Long- 
Low" observations of SNe la scheduled in G03 will be 
used to make firm identifications for the emission features 
in this band. 

3.2. SN 2005df Spectrum 

We begin our discussion of the MIR spectra with the 
second object observed, SN 2005df, as the S/N of this 
spectrum is considerably higher than that of SN 2003hv, 
and greatly aids the identification of spectral features. 
The reduced mid-infrared spectrum of SN 2005df on 14 
Dec. 2005 is shown in Fig. [51 The effective spectral 
resolution varies over the range 2500-4500 km s~^, with 
the lowest resolution occurring at the blue end of each 
order. 

Line identifications were made by searching The 
Atomic Line List^^ for strong nebular transitions at 
or near the ground state. Candidate line transitions 
were then further limited by the selection of atomic 
species with large predicted abundances in models of 
SNe la. Specific ally, we re ferred t o the 1-D deflagra - 
tion model W7 (|Nomoto. T hielema nn, fc Yokol I1984D . 
and the 1-D delayed -detonation (DD) model 5p02822.16 
(jHoflich et al.|[200^ hereafter H02) which has a moder- 
ately subluminous synthetic light curve similar to that of 
SN 1986G. Wherever possible, line identifications were 
tested by comparing the emission-line profiles of lines 
from the same atomic species. The resulting line identi- 
fications are marked in Figure[21 and are listed in Tabled 

3.2.1. Iron-Peak Elements 

The spectrum of SN 2005df is, unsurprisingly, domi- 
nated by iron-peak elements, showing strong Co emis- 
sion and also several Ni lines. Observed emission-line 
profiles for the Co and Ni lines are shown in Figures [3] 
and 21 respectively. Both of the Co lines appear to be 
fairly centered relative to the host-galaxy rest frame. The 
11.89 [Co III] line is by far the brightest feature 
in the spectrum and has a full-width at half-maximum 
(FWHM) of around 8000 km s^^ Although the S/N 
is much lower for the 10.52 /Ltm [Co II] line, the emis- 
sion profile is not a great match to that of the [Co III] 
11.89 fim line, as it appears to exhibit a broader base and, 
perhaps, a slightly narrower core. However, the [Co II] 
10.52 /im emission line is probably blended with [S IV]. 
In DD models (and W7), the S-rich ejecta are predicted 
to have much larger radial velocities than the radioac- 
tive Co-rich ejecta, making it tempting to suggest that 
the core of the 10.5 feature is largely due to [Co II] 
and the broader base due to [S IV]. 

In contrast to the Co emission, the Ni lines show a kine- 
matic offset (of order 2000 km s~^) relative to the host- 
galaxy rest frame. Unfortunately, two of the three Ni 
lines are somewhat blended with the wings of the strong 
Ar emission lines, making it difficult to compare the line 
profiles, though it appears the [Ni IV] emission line is per- 
haps slightly broader than the [Ni III] lines. The 11.00 
/im [Ni III] Hue has a FWHM around 4000-5000 km s^i. 
Finally, we note that there is no convincing detection of 
the 6.63 /im [Ni II] line, although there is some possible 
excess emission at the correct location just at the blue 

^® |http : //www ■ pa . uky . edu/$\sim$peter/at omic/ 1 



edge of the [Ar II] emission feature. The rest frame wave- 
length of this [Ni II] line is labelled in Figure^ Note that 
the 1-pixel spike near the label is clearly just noise. 

3.2.2. Intermediate-Mass Elements 

As well as the [S IV] emission likely contributing to 
the 10.52 fiui [Co II] line, the spectrum contains further 
emission from moderate- mass elements in the form of two 
Ar emission lines showing an unusual double-peaked line 
profile (Fig. [5]). As we will discuss in 14.31 such a profile 
strongly suggests a non-spherical but perhaps nearly ax- 
isymmetric distribution of Ar emission. Furthermore, as 
with the Ni emission, the Ar lines appear to be redshifted 
relative to the rest frame of the host, with the red and 
blue emission peaks appearing around -1-6300 and —4300 
km s~^, respectively. Given the dramatically different 
line profiles, it is difficult to compare the width of the 
Ar lines with those of Co and Ni. (This is further com- 
plicated by the fact that both Ar lines are blended on 
at least one side with other emission features.) However, 
taking the 8.99 /xm line as the cleaner of the pair, we 
can say that the full-width near zero-intensity (FWZI) 
is at least 23000 km s~^, compared to 16000 km for 
[Co III] 11.89 urn, and 8000 km s^^ for [Ni III] 11.00 urn. 

3.2.3. Molecular Emission? 

It is interesting to note that between the blended emis- 
sion features near 7 and 9 /zm, the flux does not drop 
away to zero as it does shortward of 6 /im, near 10 
/im, and longward of 13 fim. Instead the emission levels 
out to about 0.2 mJy in the 7.7-8.1 fim region. This 
wavelength region coincides with that of the fundamen- 
tal (Av = 1) ro- vibrational bands of silicon monox- 
ide (SiO). SiO molecular emission h as been detected in 
core-collapse supernovae (SN 1987A: lAitkeii et al.llT988l : 
Wooden 1989: 'Roc he. Aitken. fc Smithlll991l: S"N 2005af: 
Kotak et al. 2006), but molecular emission has never 
been detected in a thermonuclear supernova. However, 
the identification of SiO in SN 2005df should be con- 
sidered speculative at best, as the excess emission could 
easily be a blend of ot her unidentified faint features. 

fHoflich e t al.l (119951 ) studied the problem of molecule 
formation in detonation models of thermonuclear super- 
novae and concluded that CO might form in sublumi- 
nous SNe la, and SiO might form in SNe la with a very 
low ^^Ni yield. There is no evidence for CO fundamen- 
tal emission at the blue end of the spectrum, nor in the 
4.5 /im photometry of SN 2005df. On the other hand, 
the C/O-rich ejecta in DD models is much farther out 
than the layer containing Si and O. (In the 5p02822.16 
model, the Si and O layers overlap at a radial velocity of 
around 12500 km s~^, while C is found predominantly 
beyond 19000 km s"^) Since molecule formation rates 
are highly density dependent, it might perhaps be possi- 
ble to form SiO without forming CO. 

3.2.4. Unidentified Features and Alternative Line IDs 

A few apparently real emission features in the spec- 
trum of SN 2005df remain unidentified. Two unresolved 
features with similar brightness appear to the red of the 
strong [Co III] feature at 12.5 and 12.8 /im. The redder 
of the two features might be residual background [Ne II] 
12.8 /im emission. Indeed, examination of the raw 2-D 
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spectra does show strong [Ne II] emission, which might 
leave a small residual feature if the background is chang- 
ing near the supernova. However, we note that the 11.3 
/im [Ni I] emission is much stronger in the background, 
and there is no comparable residual emission seen in 
the supernova spectrum. It is perhaps possible that the 
residual 11.3 /im [Ni I] background could be masked by 
the [Ni III] 11.00 nm feature, though comparing the SN 
spectrum with the spectrum of nearby background emis- 
sion suggests that the two should be at least partially 
resolved. In any case, [Ne II] can only account for one of 
the two narrow features. 

Given the appearance of the Ar emission-line profile, 
it is perhaps tempting to suggest that the two features 
are, in fact, a single feature with a forked emission pro- 
file. On the other hand, the spacing of the peaks does 
not match those of the Ar lines, and indeed the entire 
feature is much narrower. In sum, no obvious ground- 
state line suggests itself for either the bluer line, or as 
an identification for a single forked feature. Similarly, 
an apparently broad emission feature appears near the 
red edge of the spectrum at 14.6 /im without an obvious 
identification. 

Finally, we note that it is possible to identify the emis- 
sion features near 7.4 and 11.1 /im as being due to [Ni I] 
7.507 and 11.308 /tm emission lines rather than [Ni III]. 
These lines are two of the four low-level fine-structure 
fines of [Ni I], and a third, at 14.81 /tm, could poten- 
tially match the 14.6 /im feature. (The fourth at 12.00 
/im would be lost in the [Co III] feature.) However, the 
implied ~ 6000 km s~^ blueshift of these three lines (see 
Fig. [6]), while self consistent, is more difficult to under- 
stand in the context of current SN la explosion models. 
It is perhaps possible that such a high-velocity blob of 
cold Ni could form in the context of a "co nfined deto- 
nation" scenario as proposed by iPlewa. Calder. fc LambI 
(j200l). 

3.3. SN 2003hv Spectrum 

The reduced spectrum of SN 2003hv is shown in Fig- 
ure [7] along with the SN 2005df spectrum, allowing a 
direct comparison of the two. The two spectra appear 
to show many of the same features, although the S/N 
of the SN 2003hv spectrum is much lower making defi- 
nite identifications of the fainter features difficult. (Error 
bars are included in Fig. [5] in order to help distinguish 
features and noise.) The [Co III] line which dominates 
the SN 2005df spectrum is still visible in SN 2003hv, al- 
though it is much weaker (both in absolute luminosity, 
and relative to other spectral features). This weaken- 
ing of the Co emission is not surprising, given that the 
SN 2003hv spectrum was obtained at an age of more than 
twice that of the SN 2005df spectrum, and indeed this 
supports the notion that most of the Co emission is due 
to radioactive ^^Co. 

Like SN 2005df, SN 2003hv also exhibits strong [Ar II] 
and [Ar III] emission features, and again they appear to 
exhibit complicated multi-peaked emission-line profiles 
(Fig. 111). In SN 2003hv it appears that the [Ar II] and 
[Ar III] line profiles may differ more significantly than 
in SN 2005df, with [Ar II] exhibiting narrow emission 
peaks at high velocities not seen in [Ar III] . The uniden- 
tified feature(s) near 12.5 /tm may also be present in the 
SN 2003hv spectrum, though there is no convincing de- 



tection of the 12.8 or 14.6 /tm features. 

It is possible that apparent excess emission near 11 /tm 
in the SN 2003hv spectrum might be [Ni III] or [Ni I] as 
well, though the S/N of this feature is very low. However, 
both of these identifications are upper-level transitions, 
and the ground-state transitions of neither [Ni III] nor 
[Ni I] are seen near 7 /zm. The red peak of the [Ar II] pro- 
file at 7.24 /im (appearing near -1-9000 km s~^ in Fig. [8]) 
is not consistent with being either [Ni III] 7.349 /zm at the 
same velocity as the 11 /im emission identified as [Ni III] 
11.002, nor [Ni I] 7.507 /im with the 11 /im emission iden- 
tified as [Ni I] 11.308 /tm. That there is no detection of 
an appropriate 7 /tm line in SN 2003hv suggests that the 
11 /im feature may be just noise. (On the other hand, the 
[Ni III] 7.349 /tm hue in SN 2005df is much weaker than 
expected compared to [Ni III] 11.002 /tm, perhaps be- 
cause one or both features are optically thick; see § 15.11 
It should also be noted that the region near 7.4 /im is 
where the edges of the three spectral orders overlap, and 
the noise in this region may be underestimated.) Weak 
features near 6.6 /im and 8.5 /tm may be [Ni II] 6.636 /im, 
and [Ni IV] 8.405 /im, respectively, though again the S/N 
of these putative features is quite low. 

4. EMISSION-LINE PROFILES 

The line profiles of the resolved and relatively well- 
isolated emission features in SN 2005df can be used to 
probe the chemical and ionization structure in the ejecta. 
To this end, we constructed model emission-line profiles 
based on relatively simple models of the kinematic dis- 
tribution of emission. The models contain no explicit 
radiation transfer or radiative emission calculations, but 
merely calculate the line profile for a given emissivity dis- 
tribution. In effect, we are using these models to infer the 
geometry of the emission distribution from the observed 
line profiles. The details of the line-profile calculations 
are given in the Appendix. 

4.1. Nickel 

With a FWZI « 8000 km s-\ the nickel hues in SN 
2005df (Fig. m are marginally resolved. Furthermore, 
two of the three features are blended with the wings of 
the strong Ar emission lines. At best it can be said that 
the 11 /tm [Ni III] line is consistent with the parabolic 
line profile of a sphere of constant-density emission. The 
observed line centers of these features suggest that the Ni 
emission is redshifted by roughly ^ 2000 km s""'^, relative 
to the rest frame of the host galaxy. 

4.2. Cohalt 

The 11.89 /im [Co III] line is also consistent with a 
simple parabolic line profile, although the central pixels 
suggest a somewhat flatter line core. Indeed the line is 
a good match to the profile of uniform hollow spherical 
emission distribution, with inner and outer radii corre- 
sponding to a minimum velocity Vmin — 2500 km s~^, 
and a maximum velocity Vmax — 8000 km s~^ (see 
Fig. [9|). Such a distribution is suggested by the 
5p02822.16 DD model, with the central hole correspond- 
ing to the electron-capture zone. The evidence for the 
central hole is relatively weak in the mid-IR spectrum 
of SN 2005df due to the low spectral resolution, but it 
is also observed in the line profile of the NIR [Co III] 
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1.55 fim feature around day 200 (see Fig. [TOl data from 
Gerardy et al., in prep). 

That the flat Hne core is not seen in the [Fe II] hne 
at this time probably indicates that gamma-rays from 
the radioactive ^^Co are still being trapped in the core 
ejecta and that the energy deposition in the ejecta is 
not yet completely local. Little Go is produced in the 
electron-capture zone, so the Go emission lines are flat- 
topped with both gamma-ray and positron dominated 
energy deposition. However, stable ^"^Fe is produced in 
the electron-capture zone and will contribute to the Fe 
line profiles with gamma-ray dominated non-local depo- 
sition. There is no evidence for a significant kinematic 
offset in the distribution of Go emission in SN 2005df. 

In contrast, the [Go III] 11.89 /im line in SN 2003hv 
is significantly blueshifted. This blueshift is also seen in 
the [Fe II] lines of the 3 90 d NIR spectrum of SN 2003hv 
(|Motohara et al.ir2006f ). Indeed, comparison of the two 
features in Figure [11] shows a generally good match be- 
tween the [Go III] 11.89 /im and [Fe II] 1.64 ^m line 
profiles, although the S/N of the mid-IR spectrum is low 
enough that the details of the [Go III] are quite uncer- 
tain. Still, the correspondence between the two features 
lends strong support to the notion that the late-time NIR 
[Fe II] emission indeed traces the distribution of radioac- 
tive ejecta as suggested by Hoflich ct al. (2004). 

4.3. Argon 

The distinctive two-pronged emission-line profile of the 
Ar lines suggests ring-shaped emission (see Fig. ITSl) , but 
other emission configurations are also possible. This 
emission profile indicates a significant lack of emission at 
low projected velocities. However, simply removing all 
the emission toward the center of spherically symmet- 
ric SN ejecta is not sufficient to produce such a profile. 
A hollow distribution will produce a flat-topped profile, 
but will not create a central trough. In fact, in the limit 
where the Sobolev approximation holds (i.e., large ve- 
locity gradients, and little continuum opacity), it is im- 
possible to create the observed argon profiles with any 
spherically symmetric homologously expanding model. 
Since spherically symmetric distributions can always be 
described as the addition of concentric spherical shells, 
the emission-line profile from any spherically symmet- 
ric distribution would always be a superposition of the 
"top-hat" emission-line profiles of spherical shells (see 
Appendix). 

Thus, a two-pronged emission profile as seen in the Ar 
lines of SN 2005df requires a significant deviation from 
spherical symmetry. These profiles can potentially be 
formed by both a relatively circular, flattened geometry 
viewed edge-on, or as a more prolate geometry viewed 
pole-on, with a large central hole in both cases. The ve- 
locities of the peaks are largely determined by the central 
hole. The peak-to-trough contrast is mostly determined 
by the extent of the deviation from spherical symme- 
try, with extreme geometries having deep troughs, and 
more spherical geometries showing shallower troughs. In- 
deed, the observed line profile of [Ar II] 6.99 ^m shows a 
much higher peak-to-trough contrast than that of [Ar III] 
8.99 /im, indicating that the asymmetry is much stronger 
in the lower ionization species. 

It is important to note that the observed argon line pro- 
files only require an asymmetric distribution of emission. 



They do not necessarily imply a non-spherical distribu- 
tion of argon in the ejecta, as the observed profiles could 
also be due to an asymmetric ionization structure or non- 
spherical energy deposition in the Ar-rich zone. Indeed, 
as discussed above, the [Ar III] emission does appear to 
be closer to spherical symmetry than the [Ar II] emission, 
and might suggest that the asymmetry is at least in part 
an ionization effect. However such and effect could also 
perhaps result from a non-spherical chemical or density 
distribution in the ejecta, and this can not be ruled out 
with the analysis presented here. 

Here we consider scenarios for both an edge-on ring-like 
symmetry and a pole-on prolate symmetry, and consider 
a potential physical mechanism for both cases. 

4.3.1. Ring Symmetry: the "Magnetic Field Model" 

As the supernova ejecta become optically thin to 
gamma-rays, the energy deposition from the decay of 
^^Go becomes dominated by positrons a nd fast electrons. 
The p ositrons are likely trapped quickly (jSoUerman et al.l 
120041 although see references therein for opposing view- 
points) but, as discussed by I Hoflich et al.l (|2004 here- 
after H04), the resulting fast electrons could travel quite 
a distance before thcrmalizing. However, in order to 
explain the flat-topped [Fe II] emission-line profiles ob- 
served in the late-time NIR spectrum of SN 2003du, the 
energy deposition from ^^Go needs to be trapped locally, 
leading H04 to suggest that a weak magnetic field was 
present in the SN ejecta. The efficiency of such magnetic 
trapping is not necessarily going to be spherically sym- 
metric and might lead to an angular dependence in the 
energy deposition from positron decay, and thus in the 
ionization structure. 

Such an effect would only be important at late epochs, 
when the ejecta become optically thin to gamma-rays 
from the ^®Co zone. In fact, the 135 d epoch of the 
SN 2005df spectrum is probably a bit early for such an 
effect, and indeed the 200 d NIR spectrum still shows 
strongly peaked [Fe II] emission (Gerardy et al. in prep; 
Fig. [TT|) . However, the Ar appears at much larger veloc- 
ities and would become optically thin more rapidly than 
the core iron. Thus, it might be possible for positron de- 
position to be important in the Ar zone while gamma-ray 
deposition still dominates in the Ni/Fe core. 

As a mock-up of such a situation we considered emis- 
sion models with a latitudinal angular profile given by 

^-(i/,) = ^"e-'^/'^" (1) 

and a power-law radial profile 

p{r) = r"^ (2) 

truncated at a maximum radius Rmax = Vmaxt- The in- 
ner edge of the Ar emission was modelled by setting the 
flux to zero within a spherical hole of radius r^, and offset 
from the center along the line of sight by Xoff- (In the 
formalism described in the Appendix, this was achieved 
by setting the ^(r, tp, 0) profile function to zero anywhere 
within this sphere.) The location of this hole is mainly 
constrained by the peaks in the line profile, which appear 
to be asymmetric in both velocity and strength. The best 
fits to these peaks were found by using a hole radius cor- 
responding to a velocity of 5000 km s^^. Just offsetting 
the sphere by 1000 km s^^ to match the observed ve- 
locity asymmetry produces too much asymmetry in the 



Mid-IR Spectra of SN 2003hv and SN 2005df 



7 



peak fluxes. Instead, we found a better match by using a 
500 km s~^ offset for the sphere, and then applying a fur- 
ther 500 km s~^ redshift to the entire line profile, perhaps 
suggesting that the Ar emission as a whole is somewhat 
redshifted, and that the central region is slightly more 
asymmetric. 

We then tried to find a set of models that would match 
both Ar line profiles with only a change in the angular 
distribution "^{ip). The parameter space is large enough 
that such "fits" are by no means unique, but a reasonable 
match (Fig. [T2)) was found for Vmax — 12000 km s~^, 
a = 1.5, (3 = 2, and Vo = 0.25° and 12° for [Ar II] and 
[Ar III] , respectively. 

The latitudinal profile functions for the two models are 
shown in Figure [131 The [Ar II] model has a very flat- 
tened distribution only covering a few degrees near the 
equator. This is driven primarily by the large peak-to- 
trough contrast, but there is some cross-talk with the 
other parameters. In particular, it is possible to fit the 
observed peak-to-trough ratio with a wider angular de- 
pendence by using a narrower radial profile (either by 
using a steeper radial power law or a lower maximum ve- 
locity) , although this tends to make it more difficult to fit 
the observed width of the wings of the Ar profiles. The 
narrowness of the angular profile apparently required 
to reproduce the observed [Ar II] line profile does cast 
some doubt as to whether such an effect could be caused 
by the asymmetric magnetic trapping described above. 
This raises the question of whether other effects (inter- 
action with an accretion disk, for example) could create 
the required extremely flattened equatorial ring-shaped 
emission. Regardless, any proposed physical mechanism 
must be able to simultaneously explain both the flat- 
tened Ar II and Ar III distributions and the essentially 
spherical Co III distribution observed in SN 2005df. 

We note that there is nothing special about the spe- 
cific form used to describe the emission here, and similar 
results can be obtained with other flattened, relatively 
circular distributions. For example, the line emission can 
also be modelled as oblate ellipsoids with a central spher- 
ical hole viewed edge-on. This leads to similar results. 
In particular, the observed peak-to-trough contrast in 
[Ar II] again requires quite a flat ellipsoidal distribution 
(major- to-minor axis ratios greater than 10:1). 

4.3.2. Prolate Symmetry: the "Off-Center DD Model" 

The argon emission-line proflles can also be reasonably 
reproduced by a prolate emission distribution with a cen- 
tral hole, and viewed more pole-on rather than edge-on 
as with the ring-shaped distribution. The model line pro- 
files shown in Figure [14] result from emission uniformly 
distributed within prolate ellipsoids with a slightly off- 
center spherical hole in the middle. 

The emission regions are shown to scale in Figure [TSl 
Both models have a major axis corresponding to an ex- 
pansion velocity of 12000 km s~^ in order to match the 
total width of the Ar features. The minor axis was then 
chosen to match the observed peak-to-trough contrast. 
The [Ar II] model thus requires a narrower (more pro- 
late) emission ellipsoid, corresponding to 6500 km s^^, 
versus 8000 km a'^ for [Ar III]. 

As with the ring proflle, the central hole is largely de- 
termined by the velocity and relative strength of the two 
peaks and again has a radius corresponding to a velocity 



of 5000 km s^^. But for these models the hole is shifted 
1000 km s~^ to the back side of the SN envelope (rather 
than 500 km s~^ for the ring models). Again, as with the 
ring proflles, the entire emission line has been redshifted 
by a further 500 km s~^. 

Such an emission model is motivated by off-center DD 
models which produce an asymmetric distribution in the 
innermost ejecta. Figure [T6| shows the argon distribu- 
tion from a 2-D hydrodynamic explosion model of an 
off-center delayed deto nation. (T h is mo del is the same 
as that described by .Fesen et al.l (|2006D , and uses the 
same color scale as their Figure 7, with red being the 
mass- fraction peak and blue being zero mass-fraction.) 

The DD model is generally more spherical than our el- 
lipsoidal geometry, but does show a pronounced crescent- 
shaped peak on one side. This emission geometry would 
only produce a single emission peak on a relatively broad 
and flat distribution. However, it is worth noting that 
the deflagration in this model was calculated in 1-D, and 
therefore strictly maintained the spherical symmetry of 
the core region prior to detonation. It is perhaps possible 
that a more realistic treatment of the deflagration might 
produce both a more asymmetric chemical distribution, 
and could perhaps produce a second argon peak, depend- 
ing on the details of how the detonation propagates past 
the deflagration core. 

5. ABUNDANCES 

The Ni, Co, and Ar features observed in the SN 2005df 
spectrum all arise from low-level fine-structure transi- 
tions. Most are ground-state transitions from upper lev- 
els with excitation temperatures of 1000-2000 K, making 
them relatively insensitive to temperature. (The one ex- 
ception is the 11.002 /im [Ni III] line, for which the lower 
level is the upper level of the 7.349 /im line, and the upper 
level has an excitation temperature near 3300 K. Conse- 
quently, this line is somewhat more sensitive to temper- 
ature.) Furthermore, the low excitation temperature of 
the levels means that the stimulated emission correction 
is a large effect, which tends to keep the optical depth of 
these features relatively low. Finally, in the early nebu- 
lar phase, the electron density in the ejecta can be well 
above the critical densities of the strong lines, which can 
make even simple LTE derived estimates of the ejecta 
abundances relatively accurate. Indeed, to some extent 
this can even be checked (see § 5.5) by using the derived 
abundances and observed line profiles to observationally 
infer the density in the ejecta. 

In practice, a fully-NLTE nebular emission model is 
probably still needed for robust abundance measure- 
ments, but such an analysis is beyond the scope of 
this paper. We hope to include such an analysis in a 
subsequent paper using the entire Spitzer dataset from 
SN 2005df as a constraint. Here we will simply use the 
observed line emission and an LTE approximation to es- 
timate the abundances. Since we are extrapolating the 
LTE abundances from ground-state transitions, NLTE 
corrections to these abundances will probably tend to 
decrease the total mass, as the upper levels deplete rela- 
tive to their LTE populations. 

Calculated abundances are given in Table [H For each 
line included in the table we have inferred an observed 
mass from the number of atoms in the upper level of 
the transition as implied by the line flux and the atomic 
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transition probabilities (values and references given in 
Table [2]) . Clearly this is an absolute lower limit to the 
abundances. Assuming LTE level populations, the total 
mass of the ionization species can be calculated using the 
population of the upper level of the transition and the 
atomic partition function. 

Note that i n the following discussion we adopt the 
iHamuvl (|200l distance of 18.7 Mpc to SN 2005df un- 
less otherwi se stated. To s c ale th e results to the 16 Mpc 
distance of 'S chmidt et all ()1994D or the 14.4 Mpc dis- 
tance of LEDA divide the masses and optical depths for 
SN 2005df by a factor of 1.4 or 1.7 respectively. 

5.1. Stable Nickel 

For Ni, we use the partition functions of lHalenka et al.l 
()2001l ). and calculate masses for temperatures of 1000, 
6000, and 10000 K. Since the excitation temperatures 
of these lines are low, the derived mass changes little for 
the higher temperatures, but is dramatically increased at 
T = 1000 K. Such a temperature is probably too low in 
the presence of local radioactive heating, but might per- 
haps be realistic if the Ni is both spatially separated from 
the radioactive Co and the Ni region is optically thin to 
gamma-rays from the radioactive zone. The former is 
certainly suggested by the line profile analysis, but the 
the NIR [Fe II] line profile (Fig. [TO]) suggests that the in- 
ner ejecta are still trapping gamma-rays. For the higher 
temperatures, the derived abundances suggest a total Ni 
mass M^i « (7 — 12) x 10~^Mq, somewhat low com- 
pared to the predictions of both W7 and the 5p02822.16 
DD model. Mat, = 6 x IO'^Mq and 7 x IO'^Mq, re- 
spectively. However, the ratio of the observed [Ni III] 
lines suggests that they are optically thick, which may 
explain the discrepancy between the derived and model 
Ni abundances. 

Still, the implied Ni masses are much larger than can 
be explained as radioactive ^^Ni at such a late epoch 
(some 15-1- half-lives after the explosion). Similarly, the 
derived Ni mass is at least an order of magnitude larger 
than could be explained by a Chandrasckhar mass of 
primordial nickel at solar abundance. This strongly sug- 
gests that the observed Ni is an electron capture product 
from high-density nuclear burning early in the explosion. 

5.2. Radioactive Cobalt 

For the Co^"*" ab undance, we use the partition function 
of lHalenkal (|1989[ ). for temperatures of 3000, 6000, and 
10000 K. In all three cases the derived Co'^^ mass is quite 
similar, Mco2+ « 7 x lO^^Mg. 

The number of ^^Co atoms at time t is given by 

^, = /2^(e-^i*-e-^^*), (3) 

where Ai and A2 are the decay constants for ^^Ni and 
^^Co, respectiv ely, and A^n is the total number of ^^Ni 
atoms at t = (|Colgate fc McKeel[l969l ). Assuming that 
Co^+ is the dominant ionization species, the 11.89 fim 
[Co III] line flux implies a relatively low ^^Ni yield, 
Afse^fj « 0.22Mq. This is well below the Msejv^ = O.6M0 
prediction of W7, and other similar estimates of the ^^Ni 
yield for normal SNe la. 

For 6000 and 10000 K, the optical depth estimates for 
the 11.89 /im line (see below) suggest that the line is rel- 



atively optically thin, and correcting for the escape prob- 
ability of 11.89 /im photons only adds another ~ 20% to 
the inferred Co^+ mass. The implied optical depth for 
3000 K is 0.8, which is non-negligible. Still, correcting 
for the escape probability a,t t — 0.8 only brings the ^®Ni 
mass up to ~ 0.32AfQ. 

It is perhaps possible that a significant mass of Co 
is in a different ionization state. However, the relative 
weakness of the 10.52 fim [Co II] line compared to the 
11.89 fim [Co III] line suggests that there is comparably 
little singly ionized Co. Unfortunately, the ground state 
[Co IV] transit ion at 15.647 iim lies just off the red edge of 
the spectrum. iRuiz-Lapuente et al.l ()1995f ) estimate that 
the Fe'^+ZFe ratio is less than 30% based on UV observa- 
tions of SNe la at an age of ~ 1 year and'SoUer man et al.l 
(j2004) similarly find that the contr i bution of Fe'^"'" should 
be low at 300 d. Finally lAjcebol (|1980D finds Fe^+ZFe 
« 0.1-0.2 at 150 days. Hence, it seems unlikely that the 
ionization could account for enough radioactive nickel to 
bring the ^^Ni yield up to match the Msejvi ~ O.6M0 
predicted by W7 and normal-luminosity DD models. 

We conclude that SN 2005df probably had a low yield 
of radioactive ^^Ni by at least a factor of two compared 
to normal SNe la. This is even more pronounced for 
the closer distance estimates for SN 2005df as both the 
^^Ni mass and the optical depth become smaller. In con- 
trast, the same calculations for SN 2003hv (a spectro- 
scopically normal SN la) yield much larger ^^Ni masses; 
Ahepfi « 0.47Mq for T = 3000 K and 6000 K, and 
Maajsii K, O.53M0 for T = 10^ K. Such masses are a 
good match to the normal-luminosity DD models of H02 
and are slightly less than the 0.6 Mq of W7, though well 
within the errors from the line flux. In fact, if we assume 
Mse^j = O.6M0, these observations would imply that 
Co^+ is by far the dominant ionization state for Co in 
SN 2003hv. 

Thus, the ^^Ni yield in SN 2005df appears low by a 
factor of ~ 2 - 4 compared to both the model expec- 
tations for normal-bright SNe la, and also to the ob- 
served ^^Ni yield from SN 2003hv, an apparently normal 
SN la. However, as we showed in § 2, it also seems that 
SN 2005df was a subluminous event, by roughly 0.5 mag 
{D = 18.7 Mpc) to 1 mag {D = UA Mpc). In fact, the 
derived ^^Ni yields for these two distances are actually 
quite close to the calculated yields for the subluminous 
DD models of H02. The 5p02822.14 model is about 1.5 
mag underluminous, and has a predicted ^^Ni yield of 
M56^j = O.15M0. This compares quite well to the re- 
sults assuming a 14.4 Mpc distance (Msejvi = O.13M0 
and 1 mag subluminous). Similarly, the 5p02822.16 
model is approximately 0.5 mag subluminous and pre- 
dicts a ^^Ni yield of A/se^j = O.27M0 which again com- 
pares relatively well to the observed results for an as- 
sumed 18.7 Mpc distance (Mse^i « O.22M0 and 0.5 
mag subluminous). Thus, it seems that while the implied 
^®Ni yield is low compared with the nominal expectations 
for normal-bright SNe la, it is actually quite a decent 
match to the predictions from subluminous DD models. 
Furthermore, if these derived ^^Ni abundances survive 
a more robust NLTE spectral analysis, and if more de- 
tailed photometric analysis confirm that SN 2005df was 
indeed subluminous, then these MIR observations pro- 
vide strong direct support for the notion that the lumi- 
nosity of SNe la is correlated with the ^^Ni yield. Note 
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that IStritzinger et al.l (|2006f ) have also recently found 
strong evidence for this correlation based on the anal- 
ysis of bolometric light-curves and detailed modelling of 
the nebular spectra of nearby SNe la. 

5.3. Argon 

For the argon abundances, we calculated atomic par- 
tition functions for Ar^ and Ar ^^ using the polynomial 
approximations of llrwini (Il981h . Again, we calculated 
LTE Ar+ and Ar^+ abundances for temperatures of 3000, 
6000, and IC* K, and the results are relatively insensi- 
tive to temperature, yielding M^r+ « 9 x 10~^Mq and 
M^r2+ w 1.5 X 1O"^M0. The resulting Ar2+/Ar+ ratio is 
thus approximately 2:1. However, the line-profile anal- 
ysis suggests that the [Ar III] emission is coming from 
a larger region. In fact, assuming the emission is dis- 
tributed as in our prolate emission model, the ionization 
is approximately evenly split in the regions with Ar+. 
The total mass in singly ionized and doubly ionized ar- 
gon, Mat ~ 2.4 X IO^^A/q, is close to the predictions 
of W7 and the 5p02822.16 DD model, which predict 
Mat « 2 X lO^^M© and 4 x lO^^Mg, respectively. 

5.4. Optical Depth 

For homolog qusly expa nding ejecta, the Sobolev opti- 
cal depth (e.g.. [ShulfT99ll) is given by 
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where ni is the number density of atoms in the lower 
level and gm gi are the statistical weights of the upper 
and lower levels, respectively. Assuming an LTE relation 
between the upper and lower level populations given by 
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the Sobolev optical depth becomes 
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From this, an LTE optical depth can be calculated (for a 
given temperature) from the total number of upper-level 
atoms implied by the line flux, converted to a number 
density by a judicious assumption about the distribution 
based on the line-profile analysis. 

The logic of such a calculation can be slightly circular, 
in that the total number of emitting atoms is inferred 
from the observed flux assuming that the line is opti- 
cally thin. This is then converted to a density to infer an 
optical depth. If the calculated optical depth is indeed 
low, then the estimate from such an analysis is proba- 
bly valid. However, if the calculated optical depth ap- 
proaches unity (or beyond), then the actual number of 
emitting atoms will be underestimated, and the actual 
optical depth will be higher, perhaps significantly so. 

For each of the abundance estimates in Table [3l we 
have also calculated an LTE optical depth. For the Ni 
lines, we assumed that the Ni was uniformly distributed 
in a sphere with a radius corresponding to an expansion 
velocity of 4000 km s^^. For the Co lines, we also as- 
sumed a uniform distribution, but this time in a hollow 
spherical volume with an inner radius corresponding to 
2500 km s~^ and an outer radius of 8000 km s~^. 



This analysis implies that the Ni lines are optically thin 
at 6000 and 10000 K, but optically thick at 1000 K. The 
observed flux ratio of the two [Ni III] lines suggests that 
these lines may indeed be optically thick. In principle, 
the [Ni III] 11.002/7.349 /im line ratio can be used as a 
temperature probe (at least for low temperatures). The 
obse rved line ratio, .6 ±0.2, is much larger than it should 
be. iBautistal (|200H ) has calculated the 11.00/7.35 /xm 
line ratio for a large range of densities and temperatures, 
and the ratio is always less than 0.2. While the observed 
line ratio is not w e ll con strained, making it consistent 
with the IBautistal (|2001h calculations would require a 
rather large change in the observed line fluxes. Indeed, 
taken at face value, the upper level populations inferred 
from the line ratio imply inverted level populations for 
the 11.002 /im line. It seems more likely that the Ni lines 
are simply optically thick. 

The optical depth of the Co lines in SN 2005df is rel- 
atively low for all the abundance estimates, except for 
T = 3000 K where it is beginning to become an impor- 
tant effect. Still, even at 3000 K, the Sobolev escape 
probability for the [Co III] 11.89 /xm photons is around 
70%, so it is unlikely that a large amount of [Co III] is 
being missed. If the temperature in the Co zone was as 
low as 1000 K, then the line would be optically thick, but 
it seems unlikely that the temperature has become so low 
in the ^^Co region 135 days after the explosion. The op- 
tical depth of the [Co III] 11.89 //m hue in SN 2003hv is 
negligible for all temperatures. 

Optical depth calculations are slightly more compli- 
cated for the Ar features, as the densities implied by the 
line emission depend on the assumptions made about 
the emission distribution. For the optical depths quoted 
in Table |3l we assumed the emission was distributed as 
in our prolate line-profile model. The resulting optical 
depths are small, suggesting that the Ar lines are indeed 
optically thin. 

5.5. Electron Density 

One check on the robustness of the derived abundances 
would be to see if the implied electron density is above 
the critical density for the observed ions. If it is, then the 
collisional processes will dominate and tend to drive the 
level populations toward LTE, and thus our LTE mass 
estimates would be relatively valid. 

Unfortunately, there are almost no published collision 
strengths for doubly ionized cobalt, so we used the ap- 
proximation 



ij — Logigj, 



(7) 



with 



0.2 (|A 2cel rod| [1985 



Graham. Wright, fc Longmord 119871 : iBowers et al.l 
1997f ). This yields a typical critical density 



ncrit ~ 10^ cm^'^ for the strong optical, NIR, 
and MIR [Co III] transiti ons. Published collision 
strengths for Ar"*" and Ar^+ (iPelan fc BerringtonI 119951 : 
iGalavis. Mendoza. fc ZeippenI I1995F give critical densi- 
ties ricrit ~ 3 X 10^ cm~^ for both of the Ar transitions 
seen in SN 2005df and SN 2003hv. 

For SN 2005df, the Co^"*" number density as calculated 
for the optical depth estimates is nco^+ ~ 1.4 x 10^ cm~'^. 
Assuming that Co^+ is the dominant ionization species, 
and an epoch of 135 days. Equation |5| implies that the 
total density of ^^Ni, ^^Co, and ^^Fe ions in the Co- 
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emitting zone is nsejvi + nseco + nsep^. « 4.5 x 10^ cm . 
In the 5p02822.16 DD model, the mass fraction of ^^^Ni in 
the 5^Ni-rich zone is 0.6-0.8, so the total ion density in 
this region should be about 30% higher than the density 
of radioactive ejecta products, Uion ~ 5.5 x 10^ cm~'^. 

The total electron density will depend on the average 
ionization state of ions in the radioactive zone. Ioniza- 
tion models for SNe la suggest that Fe+ should dom- 
inate in the ^^Co zone at an age of around 300 days 
(iRuiz-Lapuente et al.l ll995: Liu. Jefferv. fc Schultzlll997l : 
iSollerman et al.ll2004D . This would give an electron den- 
sity He ~ Uion, and thus even this relatively low ion- 
ization is enough to put the electron density above the 
critical densities for the Ar lines, and near the critical 
densities for Co. Furthermore, at this earlier epoch the 
ionization is likely to be higher, which could drive the 
electron density over the critical density. 

This suggests that the LTE approximation is not too 
bad, at least for the Co and Ar abundance in SN 2005df, 
and that NLTE level populations arc unlikely to explain 
to low ^^Ni abundance in the supernova. The same calcu- 
lation for SN 2003hv yields a significantly lower ion den- 
sity riion « 5 X lO'' cm~'^, and suggests that the NLTE 
corrections to the abundance may not be small in the 
older supernova. 

6. DISCUSSION 

One of the most controversial topics of discussion in 
SN la physics is the question of the exact nature of the 
thermonuclear burning front. In particular, does it prop- 
agate primarily as a subsonic deflagration front or as 
a near-sonic or supersonic detonation? The distinction 
is important because WDs undergoing deflagration have 
time to react to the release of nuclear energy and evolve 
dynamically during the supernova explosion, while deto- 
nations proceed through the star too rapidly for dynam- 
ical evolution to occur. 

Pure detonation models have long been ruled out as 
they burn the WD progenitor alm ost entirely to iron- 
peak elements (e.g., iNomotol |1982[ ) and cannot repro- 
duce the observed intermediate-mass spectral features 
of Si, S, Mg an d O typically seen in SN la spectra 
(jFilippenkol II997I and references therein). In contrast, 
1-D pure deflagration models such as W7 result in a 
spectrum of burning products ranging from iron-peak, 
to intermediate-mass elements, and also leave a sig- 
nificant portion of the C/0 pr ogenitor imburned (e.g., 
iNomoto. Thielemann, fc Yokoil[T984) . 

A similar spectrum of burning products results 
from the d elayed-detonation scenario (i Khokhlovl II991I : 
lYamaoka~e t al. 1992), where the burning begins as a de- 
flagration but quickly transitions to a detonation. The 
deflagration serves to pre-expand the WD, lowering the 
density of the outer regions of the star so that the result- 
ing detonation produces a spectrum of burning products 
and not simply Fe-peak elements. Since the detonation 
propagates too fast for the formation of convective tur- 
bulence, the WD undergoes essentially no hydrodynamic 
evolution during that stage of the burning, and thus 
a radial chemical st ructure is produced even in mult i- 
dimensional models (jCamezo. Khokhlov. fc Oranl[2005h . 

While 1-D deflagration models produce a radi- 
ally stratified structure, 3-D deflagration models have 
shown that the burning is subject to instabilities 



which result in the formation of large-scale convec- 
tive plumes and macroscopi c mixin g of the ejecta (e.g ., 
iGamezo. Khokhlov. fc Oralll [200I iRopke et al.l [200I . 
As a result, 3-D deflagration models predict a much more 
chemically homogeneous ejecta structure. 

In the two SNe observed here, we find strong evidence 
for a layered chemical structure like that seen in DD 
models (or 1-D deflagrations like W7). The line profiles 
in SN 2005df suggest a clear radial progression from Ni 
{v < 4000 km s-i), to Co (2500 < v < 8000 km s'^) 
to Ar (5000 <v< 12000 km s^^). The detection of Ni 
in SN 2003hv is less clear, but certainly the Ar lines are 
again much broader than the observed Co line in this 
event. Thus the MIR spectra come down quite clearly 
against the sort of large-scale mixing predicted by the 
3-D defiagration models. 

Indeed, the observed radial extent of these layers in 
SN 2005df matches quite well to the predicted velocities 
from the 5p02822.16 DD model (see Fig. [HI). In this 
model, the Ar runs from ^ 5000 to 12000 km s^^, the 
radioactive ^^Ni zone from ^ 3000 to 8000 km s^^, and 
the stable Ni in the inner ~ 3000 km s~^, with a lower 
abundance "tail" out to - 4000 km s-i(H02). 

However, this innermost radial structure poses a bit of 
a problem even for DD models. The high-density elec- 
tron capture burning in the II02 DD models all takes 
place during the defiagration phase, and thus should be 
subject to the same instabilities and mixing suggested 
by the 3-D deflagration models. (The fact that the II02 
DD models do not show this is due to the fact that the 
deflagrations are calculated in 1-D rather than 3-D.) In 
the 5p02822.16 model the detonation begins at a velocity 
of 4000-5000 km s~^ (marked as a vertical dashed line 
in Fig. [T7|) . so the layers below this velocity should show 
signs of mixing. In SN 2005df, the Co does not appear to 
extend much below ^ 2500 km s""'^, and the Ni is conflned 
to a bubble of about ~ 4000 km s~^ in radius. In fact, 
the only thing that suggests any sort of mixing in the in- 
ner regions is the offset of the Ni from zero velocity, but 
this seems less like t he highly mixed structures of 3-D 
deflagration mode ls (jCamezo. Khokhlov. fc OranI 120041 : 
iRopke et al.ll2006l ). and (qualitatively at least) more like 
the buoy antly rising bubble of the "con flned detonation" 
scenario ewa. Calder. fc La^[200l . 

This is not the first set of observations that has failed 
to find evidence for this mixing in the innermost ejecta 
of SNe la, although it is probably the clearest and most 
direct constraint. Spherical DD models could reproduce 
the NIR spectral evolution of the subluminous SN 1999by 
quite well. However, mixing the core of these models to 
simulate mixing during the defiagration destroyed the 
good match to the observations (H02). Support for a 
chemically stratified ejecta structure was also found by 
iBadenes et al.l ()2006[ ). in their analysis of the X-ray spec- 
tra of the Tycho supernova remnant (SN 1572). Again, 
they found that the X-ray spectrum of Tycho could be 
quite well reproduced by modelling the interaction of the 
ejecta from a 1-D delayed detonation model with the sur- 
rounding interstellar medium. 

Late-time NIR observations of some SNe la show 
fiat-topped boxy [Fe II] line profiles, indicative of a 
central ho le in the distribution o f radioactive ejecta 
(iHofiich et~al.ll2004tlMotohara et an[2006,) . Such a radia- 
tion hole can be explained as the result of nuclear burning 
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at very high densities during the deflagration, which pro- 
duces stable Fe and Ni isotopes via electron capture (EC) 
rather than radioactive ^^Ni. Again, however, if the cen- 
tral region is mixed, then the non-radioactive electron- 
capture products will be mixed with radioactive ejecta 
and a central radiation hole would not be observed. 

Possible alternative explanations for this observed cen- 
tral hole in the radioactive ejecta remain unconvincing. 
Some 3-D deflagration models (e.g. [Ropkc 2005) predict 
a large fraction of unburned material at low velocities. If 
the chemical differentiation is strong enough, this might 
produce a central hole in the radioactive ejecta and thus 
the boxy line profiles. However, such a model would not 
explain the large mass of low-velocity stable Ni seen in 
SN 2005df. Furthermor e, mo dels of the nebular spec- 
tra based on the Ropkg l|2005[ ) explosion models predict 
strong carbon and oxygen emiss ion which is not se en in 
the late-time spectra of SNe la (|Kozma etall [20051) . In- 
teraction with the secondary star in a single-degenerate 
scenario m ight produce a nearly conical hole in the ejecta 
(|Marietta et al. 2000; Kascn ct al. 2004), but this wiU ef- 
fect all layers of the ejecta equally, or nearly so, and will 
not produce a central hole in only in the innermost layers. 
(However it is possible that this interaction might could 
perhaps explain some of the bulk kinematic offsets seen 
in the late-time NIR spectra of some SNe la as reported 
bv lHoflich erall 120041 : iMotohara et al.ll2006l ) 

The MIR observations support the EC radiation hole 
explanation for the boxy [Fe II] line proflles. In 
SN 2005df, we see direct evidence that a significant 
amount of stable nickel is formed in the innermost ejecta, 
which is confined to small radial velocities, and that there 
does appear to be a hole in the center of the distribution 
of radioactive ejecta. (These do not quite coincide in ve- 
locity space since the Co distribution is centered, while 
the Ni is slightly redshifted.) Furthermore, the [Co III] 
in SN 2003hv is significantly shifted to the blue, and 
matches the observed blueshift of the NIR [Fe II] line 
(|Motohara et al.ll2006D . which clearly supports the no- 
tion that the late-time NIR [Fe II] lines can be used to 
infer the distribution of radioactive ejecta. 

This observed lack of mixing has now been estab- 
hshed in at least four SNe la (SN 1999by, SN 2003du, 
SN 2003hv, and SN 2005df), and as a result of three en- 
tirely independent lines of inquiry (photospheric epoch 
NIR spectra, late-time [Fe II] line profiles, and mid-IR 
line profiles.) Furthermore, this constraint docs not sim- 
ply come down on the side of detonations rather than 
deflagrations. We are observing a lack of mixing in the in- 
nermost ejecta, which burns as a deflagration in both the 
DD and pure deflagration scenarios. This seems to sug- 
gest that something important is still missing from our 
understanding of SN la physics. Maybe something other 
than a deflagration is responsible for the pre-expansion 
of the WD progenitor before a detonation. Alternatively, 
perhaps there is some key piece of physics missing from 
the 3-D deflagration models which maintains the radial 
structure of the 1-D models. 

7. SUMMARY 

We have presented the first two mid-IR spectra of 
SNe la. Although such observations are currently still 
quite difficult, requiring both nearby SNe and space- 
based instruments, the mid-IR is certainly a rich source 



of important data on these events. 

The observed line profiles of SN 2005df show clear 
evidence of a highly stratified chemical structure, with 
a small Ni-rich zone near the center, and progressively 
larger surrounding ^^Co and Ar layers. That the in- 
nermost regions of the ejecta are so stratified seems to 
strongly support a detonation scenario. Such a layered 
structure is incompatible with the highly mixed chemical 
structures predicted by current 3-D deflagration mod- 
els. Furthermore, this layered chemical structure per- 
sists even in the innermost ejecta, which should con- 
tain the products of deflagration burning in both the 
DD and pure-deflagration scenario, and suggests that a 
key piece of physics may still be missing from the multi- 
dimensional models. 

The line profiles of both SN 2003hv and SN 2005df 
also show signs of deviations from spherical symme- 
try The [Co III] line in SN 2003hv was blueshifted by 
~ 3000 km s~Y cons istent with the NIR [Fe II] fines 
(jMotohara et al.l 120061 ). Such a kinematic offset in the 
radioactive ejecta can perhaps be explained as the result 
of an off-center delayed-detonation. 

In contrast, the [Co III] line in SN 2005df shows no 
significant kinematic offset, although there is an offset 
in the distribution of stable nickel in SN 2005df. This 
cannot be explained with an off-center detonation, as 
the stable Ni is formed via electron capture during the 
deflagration phase. It is perhaps somewhat reminiscent 
of the buoyantly rising bubble in the "confined deflagra- 
tion" scenario. 

The observed Ar line profiles in SN 2005df require 
an emission distribution that deviates significantly from 
spherical symmetry, especially for [Ar II]. The implied 
asymmetry is much more pronounced than that seen in 
the calcium and iron distribution in SN 1885 (Fese n et all 
[2006) ■ and is much larger than the ^ 20% asphericity im- 
plied by the strong continuum polarization in SN 1999by 
(jHoweU et alll2001h . The difference between the [Ar II] 
and [Ar III] line profiles suggests that much of this may 
be due to an asymmetric ionization structure rather than 
an asymmetric chemical structure. 

Our prolate Ar emission model for these lines is some- 
what similar to the predictions of off-center DD mod- 
els, although currently these only produce a single Ar 
peak. Currently these models also still assume a spher- 
ically symmetric deflagration, and therefore are also in- 
compatible with the observed Ni offset in SN 2005df. It 
remains to be seen whether relaxing this symmetry in an 
off-center DD model would produce a second Ar peak. 

On the other hand, there is little sign of asymmetry in 
the radioactive ^^Co in SN 2005df, so it may be difficult 
to explain all three features in an off-center DD scenario. 
Indeed, we have also shown that the Ar lines can be 
reproduced in a more ring-shaped geometry, although it 
requires a very flat distribution of [Ar II] emission. 

Finally, we have used the observed [Co III] emission 
to infer ^^Ni yields for both SNe using an LTE approxi- 
mation, and assuming that the dominant cobalt species 
is Co2+. For SN 2003hv this yields an estimate of 
Msejvi ~ O.5M0, although the derived electron density 
in the Co zone may be well below the critical density for 
this SN. In contrast, for SN 2005df, we flnd a low radioac- 
tive ^^Ni mass, Mse^j « O.13-O.22M0, and in this case 
the electron density should be near or above the critical 
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density for Co^+ emission, so the LTE estimate should be 
reasonably accurate. Furthermore, it seems unlikely that 
other ionization species could add enough to bring this 
up to the nominal Msejv, ~ O.5-O.6M0 level. SN 2005df 
appears to have been subluminous, and so this low ^^Ni 
yield seems to be direct evidence to support the notion 
that the luminosity of SNe la is primarily determined by 
the ^""Ni mass. 

This work is based on observations made with the 
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TABLE 1 
SN 2005DF MIR Photometry 



JD Obs. 


Epoch^ 




F„ (mjy) 






(245000+) 


(days) 


3.6 ^ia° 4.5 ^m'^ 


5.8 fj.m'^ 


8.0 /.tm° 16 /.tm' 


22 fim^ 


3676.5 


93 


0.36 ±0.03 0.08 ±0.03 


0.11 ± 0.07 


0.78 ±0.22 




3689.1 


105 






1.18 ±0.52 




3763.1 


179 






0.95 ± 0.54 




3774.5 


191 


0.10 ±0.03 0.03 ±0.02 


< 0.07 


0.68 ±0.17 




3952.6 


369 






< 0.31 


< 0.92 


3955.6 


372 


< 0.03 < 0.02 


< 0.11 


0.40 ±0.16 




Assuming explosion 


occurred on 


2 Aug. 2006*^ IRAC Ch. 1" IRAC Ch. 


2^* IRAC Ch. S'* 


IRAC Ch. 4' PUI "Blue"s PUI 


"Red" 






TABLE 2 












SN 2005DF Emission Lines 








Species 


Transition 


Arest(/^rn) 


Afc, (s-i) 


Line Flux'' 


6.6 


[Ni II] 


3d9 2£,5/2-3rf^ '^3/2 


6.636 


5.5 X 10^2 


< 5'^ 


6.91, 7.16 


[Ar II] 


3p5 2p°/2-3p5 2p°/^ 


6.985 


5.28 X 10~2 d 


28 ± 3 


7.40 


[Ni III] 


3d« ^F4-3d^ ^Fa 


7.349 


6.5 X 10"2 


6 ± 2'' 


7.7-8.1 


SiO 


Ay = 1 


~ 7.5-8.0 




~ 4 


8.49 


[Ni IV] 


W 4F9/2-3<i^ 4^7/2 


8.405 


5.67 X 10~2 " 


7 ±2'' 


8.90, 9.20 


[Ar III] 


3^* 3p2_3p4 3pj 


8.991 


3.1 X 10-2 f 


23 ± 2 


10.60 


[S IV] 


3p 2P°^2-3P ^^3°/2 


10.510 


7.7 X 10-3 ^ 


11 ± 2 




[Co II] 




10.521 


2.21 X 10-2 g 




11.14 


[Ni III] 


3d» 3P3-3d8 3^2 


11.002 


2.7 X 10-2 d 


3.5 ±0.5 


11.99 


[Co III] 


3d'^ a'^Fg/2-3d'^ a*i^7/2 


11.888 


2.01 X 10-2 d 


31 ± 2 


12.5 


? 






~ 2 


12.8 


? 








~ 2 


14.6 


7 








~ 2 



In units of 10 erg cm 2 s ^.'^ Flux uncertain: line is blended with stronger feature."^ IQuinet fc Le Dourneufl 
(1996).'^ Nussbaumer & Storey (1988). <= Mclcndez fc Bautis tj ||2005|) .^ NIST Atomic Spectra Database (V3.0.3) 
http : //physics .nist .gov/PhysRefData/ASD/ index .html . ' .Quinet, (il998l) . 
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TABLE 3 

Derived Abundances from Mid-IR. Emission Lines 



Line ID 


Obs. Mass (Mq) 




LTE Inferred Mass^(M0) and Ti,^ 








1000 K 


6000 K 


10000 K 



SN 2005df= 



[Ni II] 


6.636 /^m 


< 6.1 X 10" 


-4 


< 8.7 X 10-3 


< 1.6 


< 2.7 X 10-3 


< 0.1 


< 4.1 X 10-3 


< 0.06 


[Ni 111] 


7.439 fira 


7.0 X 10- 


-4 


7.1 X 10-3 


2.4 


2.4 X 10-3 


0.1 


2.7 X 10-3 


0.08 


[Ni 111] 


11.002 


1.5 X 10" 


-3 


7.8 X 10-2 


3.0 


8.8 X 10-3 


0.2 


9.0 X 10-3 


0.2 


[Ni IV] 


8.406 ^lra 


1.1 X 10" 


-3 


8.7 X 10-3 


3.6 


3.9 X 10-3 


0.2 


4.3 X 10-3 


0.1 










3000 K 




6000 K 




10000 K 




[Co 111] 


11.888 Atm 


1.9 X 10" 




6.8 X 10-^ 


0.8 


6.6 X 10-^ 


0.4 


7.5 X 10-2 


0.2 


[Ar II] 


6.985 Aim 


2.4 X 10" 


-3 


1.2 X 10-2 


0.2 


9.0 X 10-3 


0.1 


8.2 X 10-3 


0.06 


[Ar 111] 


8.991 ^lra 


4.3 X 10" 


'3 


1.7 X 10-2 


0.2 


1.5 X 10-2 


0.1 


1.5 X 10-2 


0.06 



SN 200311V 

[Colli] 11.888 Atm 4.5 x 10-3 1.7 x 10-2 0.03 1.7 x 10-2 0.01 1.8 x 10-2 0.007 

Using atomic partition functions of Hal enka et al.l l|2001l ) for Ni; IHalenlsai {1989) for Co: Hrwiiil l| 19811 ') for Ar.'' Assuming Vmax=4:000 km 
s-^ for Ni; timi„=2500 km g- ^ and t)maa;=8000 km s— ^ for Co. For Ar, the optical depths assume the density distribution of the prolate 
line profile models.'^ SN 2005df abundances and optical depths calculated with an adopted distance of 18.7 Mpc. 
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SN aOOSdf MIR Photometry 
AIRAC 2453676.5 (93 d) 
APUI 2453689 1 {105 d) 
□ /R4C 2453774.5 (191 d) 
MPUI 2453763.1 (179 d) 
QIRAC 2453955.6 (372 d) 
♦ PU/ 2453952.6 (369 d) 



1 / 






1 1 


1 E 

1 i 


] 1 




) 1 



0,6 



Log Wavelength (/^m) 



Fig. 1.— Observed MIR broadband photometry of SN 2005df at three epochs in Nov. 2005, Feb. 2006, and Aug. 2006. The spectral 
energy distribution shows little evidence of evolution except for the rapid fading of the IRAC Ch. 1 (3.6 /im) flux between the first two 
epochs. 
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2 - 
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SN 2005df 

14 Dec 2005 (t~135 d) 



[Ar II] 





SiO? 



<k\ ijllby/^^^ 



[Co III] 
I 

/I 




8 9 10 11 12 13 14 15 

Wavelength {/Jtm) 



Fig. 2. — The observed mid-infrared spectrum of SN 2005df. Wavelengths are shown as vacuum coordinates in the observer's frame, and 
are plotted on a logarithmic scale so that the observed line width is proportional to the velocity line width throughout the large wavelength 
span. See text for discussion of line identifications. 
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^ 1 

0.8 
0.6 
0.4 
0.2 



[Co III] 11.89 fxm 



I I I I I I I I I I I I I I I I I I I 
[Co II] 10.5a yum 
■+[5 IV] 10.51 Ajm 



-20 -10 10 

V, (10^ km s"') 



20 



Fig. 3. — The observed line emission from Co lines in the mid-IR spectrum of SN 2005df on day ~ 135, plotted in velocity space relative 
to the rest frame of the host galaxy. 




-10 10 

V, (10^ km s"') 

Fig. 4. — The observed line emission from Ni lines in the mid-IR spectrum of SN 2005df on day ~ 135, plotted in velocity space relative 
to the rest frame of the host galaxy. 
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0.5 



[Ar II] 6.99 




-10 10 

V, (10^ km s"') 



Fig. 5. — The observed line emission from Ar lines in the mid-IR spectrum of SN 2005df on day ~ 135, plotted in velocity space relative 
to the rest frame of the host galaxy. The central vertical dotted line indicates emission at rest relative to the host, while the outer lines 
indicate the approximate positions of the peaJcs in the Ar emission at —4300 and +6300 km s~^. 



.1 

■ 


[J 

Ni I] 


[\^ 

n 

7.507 






1 1 1 1 1 1 1 1 


- 


-J' 
Ni I] 


11.307 


-h- 

n 


r- 

m 


1 

r 




i 


I] 1 


4.814 


-1— 


n 





-20 -10 10 20 

V, (10^ km s"') 



Fig. 6. — Observed velocity profiles of three emission features in the mid-ER spectrum of SN 2005df on day ~ 135 that could potentially 
be consistent with [Ni I] emission at high velocity. 
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1 - 



SN 2005df (+1.2) 

14 Dec 2005 (t~135 d) 

SN 2003hv (x2) 

1 Sep 2004 (t~375 d) 

[Ar II] 



[Co III] 



I 



> 

+ 




8 9 10 

Wavelength (jim) 



Fig. 7. — The observed mid-IR spectrum of SN 2003hv compared with that of SN 2005df. Wavelengths are shown as vacuum coordinates 
in the observer's frame, and are plotted on a logarithmic scale so that the observed line width is proportional to the velocity line width 
throughout the large wavelength sijaii. Error bars (based on the output error spectra from the individual order extractions, propagated 
through the weighted averaging of the orders) are included to help distinguish noise from features in the low S/N spectrum of SN 2003hv. 
See text for discussion of line identifications. 
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Fig. 8. — The observed line emission from Ar lines in the mid-IR spectrum SN 20031iv on day ~ 375, plotted in velocity space relative 
to the rest frame of the host galaxy. The central vertical dotted line indicates emission at rest relative to the host. 




V, (10° km s ') 



Fig. 9. — The observed [Co III] 11.89 nm emission-line profile for SN 2005df on day ~ 135, compared to model line emission from a 
hollow uniform spherical distribution of emission with an inner radius corresponding to a minimum velocity Vmin — 

2500 km s-i, and an 

outer radius corresponding to a maximum velocity Vumx 

= 8000 km s-^ 
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-10 10 

(10^ km s"') 

Fig. 10. — The observed [Co III] 11.89 fim emission-line profile from SN 2005df on day 135, compared with the near-IR emission-line 
profiles of [Co III] 1.55 fim and [Fe II] 1.64 fira near day 200 (data from Gerardy et al., in prep.) 
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Fig. 11. — The observed [Co III] 11.89 fira emission-line profi le from SN 2003hv at 375 d, compared with the 410 d near-IR emission-line 
profile of [Fe II] 1.64 ^tm (NIR data from lMotohara et al.ll200g '). 
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Fig. 12. — The observed Ar line profiles in SN 2005df on day 135 compared with calculated emission-line pr ofile s for an edge-on 
equatorial ring model. The emission models for the two profiles differ only in the latitudinal profiles used (see Figure [TSj, with the [Ar II] 
emission confined to a narrow region near the equator and the [Ar III] emission coming from a much larger range of latitudes. (See text 
for further details.) 




V' (degrees) 

Fig. 13. — The latitudinal angular profile functions used for the models shown in Figure [12] 



22 



Gerardy et al. 




V, (lo'' km s"') 

Fig. 14. — The observed Ar line profiles in SN 2005df on day ^ 135 compared with calculated emission-line profiles for a pole-on prolate 
emission geometry with an off-center spherical hole near the middle. (See text for details.) 




Fig. 15. — Schematic of the model emission used to create the line profiles in Figure [141 Emission is uniformly distributed in prolate 
ellipsoids with a central spherical hole offset along the line of sight. The figure is to scale with the geometry used for the Ar line profiles 
in Figure [141 
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Fig. 16. — Argon dist r ibutio n from a 2-D model for an off-center delayed-detonation model. The model shown is the same as that 
described by IFesen et al.l I I2006I) . and the color scale is the same as their Figure 7, with red corresponding to the mass-fraction peak, 
and blue corresponding to zero mass-fraction. The box dimensions correspond to velocities of ±27500 km s~^. The model shows a 
"crescent-shaped" density peak which is perhaps qualitatively similar to one lobe of our prolate Ar emission model (see text for discussion.) 
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Fig. 17. — Abundances for the 5p02822.16 delayed detonation model (H02) plotted versus expansion velocity. For simplicity, many 
of the less abundant species have not been plotted. The vertical dashed line represents the approximate location of the transition from 
deflagration to detonation. Material to the left of this line is formed during the deflagration phase, while material to the right is the product 
of detonation. Note that the small zone of Si/S-rich material near 5000 km s~^ is an artifact of the imposed spherical symmetry of this 
model and is not seen in multi-dimensional DD models. 
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APPENDIX 

CALCULATION OF NEBULAR EMISSION-LINE PROFILES 

For nebular supernova emission with negligible sources of continuum opacity and background emission, line profiles 
can be calculated in the Sobolev approximation by adding up the contributions from differential volume elements in 
the SN envelope. Adopting the notation of Shu (1991), the contribution to the luminosity emitted toward the line of 
sight k from material with specific intensity I^, moving at (large) velocity u is given by 

dL,(k)=4^/, (^)k.V(k-u) (^^^ . (Al) 

For a homologously expanding envelope with u = -jf, we have k • V(k • u) = t^^ . Assuming that the line emission is 
intrinsically narrow, and placing the line of sight along the x axis, then line emission with a rest-frame frequency fo 
will only contribute to the line profile at observed frequency v \i v ^ v^]^! + ^] = ^^o [l + ^] • Thus 

dV\ ( dx\ ( dV\ 



J \dv J \dx J 

In other words, the luminosity of the line profile at a given observed radial velocity Vx = x/t can be calculated by 
simply integrating the specific intensity of the emitting material in planar volume elements perpendicular to the line 
of sight. Thus, for a given spatial distribution of I,^, the line profile is, in principle, easy to calculate. 

In practice, however, supernovae are easier to describe as close to spherical or axisymmetric, so it is useful to convert 
this integration over a plane perpendicular to the line of sight, to an integration over spherical coordinates centered 
on the supernova and inclined at an angle i to the line of sight (such that i = implies viewing the supernova "pole- 
on" and i = 90° corresponds to viewing the supernova equatorially or "edge-on"). In general, the specific intensity 
distribution can be decomposed into the product of three components, 

/.(r, V, 0) = p(r, V, 0)*(V', 0)'i>(0), (A3) 

where r and 4> are the standard radial and longitudinal coordinates and ip is simply the complement of the standard 
latitudinal coordinate 9. For highly symmetrical distributions, the distribution will decouple further with ^'(V', 0) = 
'^{ip) for axisymmetric distributions, and p{r,^,(f)) = p{r) for spherical symmetry. (Note: we chose to use ip rather 
than 9 because it slightly simplifies the case of equatorial viewing. This derivation can be relatively easily converted 
to the standard spherical coordinates by a judicious exchange of sines and cosines and swapping integration limits.) 

It is useful to think of the differential volume elements as circular rings of constant r and tp. Each ring will have a 
radius R = rcostp, inclination i, and the x coordinate of the center will be at xq — rsin-0cosi. For a given radial 
velocity Vr — x/t, the ring will only contribute the line profile at (at most) two points = ztcj)' (except for the special 



case of i — 0), where 0' = cos ^ 
Equations IA2I and IA31 we have 



r cos sin i 



The differential volume element dV = r cos ipdrdipd(p, and so combining 



dL^{x)^A'Kp{r,iP,(t))^{'4},(t>)<^{4>)[5{4> - 4>') + 5{(P + (P')y cosipdrd^ -- 

\dx 

COS zbdvd'ijj 

= 47rp(r, V^, 0)vE'(V, <^)$(0)[(5(0 - 0') + <5(0 + 0')]- 



cos^ ip sin^ i — (x — xq)'^ 

The total contribution to the line flux at a given radial velocity can thus be found by integrating over those rings that 
intersect the x = v/t plane. Hence, 

L,{x^Vr/t)^ r^p{r) ^(r, V>)$*(x, V^, r) (A4) 

■Ix Ji/;„i„(r,a;) \/ r'^ cos^ tp sin^ i ~ {x — xo{r, ip)y 



$*(a;,V,r) = [$(</>' (x, V, 0) + $(-(^'(x, V, r))] 



$*(a;,V,r) = / $((/))d0 



2-K 



where 

for i ^ and 



for 1 = 0. The latitudinal limits of integration are more complicated, and for a given r depend on the value of x relative 
to xn = rcosi and xg = — rcosi, the x-coordinates of the tp = ±90° poles on a sphere of radius r. For a given value 
of r, the latitudinal integration limits are 

For xn < X < r: 

V^mm(^; ■^) — 2' ^ ^ 

'4'max{r,x) ^ ^ + a-i 
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Fig. 18. — Calculated emission-line profiles for three simple emission geometries: a uniform sphere (left), a uniform thin spherical shell 
(center), and a uniform thin circular ring (right). 



For < X < Xn- 



For xs < X < 0: 



For — r < X < Xs'. 



^min ij'^ X^ 
^max X^ 



where a — cos^"'^(^). 

For the hne profiles presented in this work, the radial, latitudinal, and longitudinal specific intensity distribution 
functions (Eqn IA3|) were coded in Python and Equation IA4I numerically integrated using the quadpack algorithms 
in the SciPy^^ software package. The resulting line profiles from simple emission geometries are shown in Figure [TSl 
A homogeneous sphere of emission produces the well-known parabolic line-profile shape, while a thin, uniformly 
emitting spherical shell produces a "top-hat" line profile, and a thin, uniformly emitting ring produces a forked profile 
[i^ oc (1 — {v / VmaxY)~^^'^]- Sincc the volume integration in Equation IA4I is "additive," it is often useful to think of 
the emission from more complicated geometries as linear combinations of simpler emission geometries. Thus a thick 
sphere of emission with a central hole can be built up either by adding "top-hat" line profiles from progressively larger 
shells, or by subtracting the line profile from a small homogeneous sphere from that of a larger homogeneous sphere. 
In both cases one ends up with a line profile that is flat (shell- like) in the core, and parabolic (sphere- like) in the wings. 



